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Abstract 



We analyzed the ^"^Nb on ""'^Mg reaction at 30 MeV/nucleon in the aim of disen- 
tanghng binary sequential decay and multifragmentation decay close to the energy 
threshold, i.e. ~ 3 MeV/nucleon. Using the backtracing technique applied to the 
statistical models GEMINI and SMM we reconstruct simulated charge, mass and ex- 
citation energy distributions and compare them to the experimental ones. We show 
that data are better described by SMM than by GEMINI in agreement with the 
fact that multifragmentation is responsible for fragment production at excitation 
energies around 3 MeV/nucleon. 

Key words: Backtracing technique, multifragmentation, binary sequential, 
excitation energy, SMM, GEMINI 
PACS: 25.70.-Z 25.70.Pq 24.10.Pa 



1 Introduction 

Heavy ion collisions at intermediate energy (20-100 MeV/nucleon) feature 
several mechanisms, such as fusion, deep inelastic collisions (DIG) and direct 
reactions present at low energy, or incomplete fusion and multifragmentation, 
characteristic of higher energies [1] . Through heavy- ion collisions, very excited 
nuclei can be produced offering the possibility to study nuclei far from their 
ground state conditions. The investigation of the thermodynamical properties 
of nuclear matter, such as temperature, density, excitation energy, can then be 
carried out. There is intense theoretical and experimental research surrounding 
the debate on the presence of a liquid-gas phase transition at the origin of the 
multifragmentation process [2]-[S]- 

At low excitation energies, the decay of hot nuclei follows the predictions of 
the statistical model for a compound nucleus, and the deexcitation occurs by 
binary sequential decays [7J,[8J. For systems excited at higher energy around 3 
MeV/nucleon |9], the multifragmentation channel opens which is interpreted 
as a manifestation of a liquid-gas phase transition (TU]. Many experimental 
clues have been collected recently [3], [11] on multifragmentation and its rela- 
tion to phase transition. Nonetheless the question is still debated. Disentan- 
gling binary sequential decay from multifragmentation close to the threshold 
should help to better define the thermodynamical properties mentioned above. 

In this aim we will study the decay of hot nuclei formed in the ^^Nb -|-"°* 
Mg reaction at 30 MeV/nucleon |T2]. Using the backtracing technique [13J- 
[15] the data will be compared with the predictions of two codes : GEMINI 
[TB] for binary sequential decay, and the Statistical Multifragmentation Model 
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(SMM) [T7] for multifragmentation. The reaction dynamics was simulated 
by Boltzmann-Nordheim-Vlasov (BNV) transport model simulations which 
describe the evolution of the one-body distribution function according to the 
nuclear mean field and including the effects of two-body collisions with the 
test-particle method |18], [19]. We used a soft equation of state {K = 200 MeV) 
and a free nucleon-nucleon collision cross-section a^N- The time evolution of 
the density contour in the reaction plane for various impact parameters was 
simulated [20j. Those simulations [12] predict the formation of a single source 
by incomplete fusion for small impact parameters [b < 4 fm). For larger impact 
parameters the mechanism showed a binary character with formation of two 
sources, quasi-projectile (QP) and quasi-target (QT). In this paper we will 
study the decay of incomplete fusion and QP sources. 



2 The Experiment 

Collisions of ^^Nb ions at 30 MeV/nucleon with a 2mg/cni? thick '^"■'^Mg target 
were studied at the GANIL facility using the INDRA multidetector. INDRA is 
constituted by 324 independent telescopes on 16 rings : the telescopes covering 
the polar angular range from 3° to 45° are comprised of an ionization chamber, 
a 300 /im thick silicon detector and a CsI(Tl) scintillator. Those covering the 
range from 45° to 176° are comprised of two layers : an ionization chamber 
and a CsI(Tl) crystal EI]-[21]. 

The reaction ^^Nb Mg was measured during a campaign of experiments 
for which INDRA was coupled to the CHIMERA ([25] - f27]) first ring, which 
covered the angular range from 1° to 3°, and some INDRA modules in the hor- 
izontal plane were replaced by silicon strip detectors. However, for the present 
reaction study, none of these additional detectors were used. INDRA detects 
charged particles and fragments with an efficiency close to 90% of the whole 
solid angle and a high granularity in order to reduce double counting (down 
to 5%). The identification thresholds are low (~ 1 MeV/nucleon). INDRA can 
measure ion charge and energy in a wide range and it can resolve masses up 
to Z=4. 

The energy calibration of the silicon detectors and ionization chambers for 
heavy fragments was obtained by elastic scattering of Ar, Ni and Xe beams on 
a gold target at incident energies from 7 to 9 MeV/nucleon. For the scintillators 
the energy calibration was accomplished with secondary beams of hydrogen 
and helium isotopes at different energies. The charge identification was realized 
by means of a seven parameter fit of the AE-E matrices [28] which well 
reproduces the form of the lines for each atomic number, Z. Unit charge 
resolution was obtained for all nuclei produced in this reaction. 
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The data for the ^^Nb Mg at 30 MeV/nucleon were recorded with an 
acquisition trigger requiring at least 5 fired telescopes. The calculated reac- 
tion cross section is cr = 3.4 barns [29j. The total measured cross-section, 
calculated from the target thickness, integrated beam current, and total num- 
ber of recorded events, corrected for dead time, is cr = 2. barns. The lower 
experimental value is due to the lack of detection at forward angles (< 3°) 
and, above all, to the acquisition trigger condition which eliminates the most 
peripheral collisions. For this reaction, the available center of mass energy is 
EcM = 572 MeV and the projectile velocity is Vpj = 7.61 cm/ns. The grazing 
angle is Ograz = 1-44°. 

3 Data Analysis 

3.1 Event Selection 
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Fig. 1. Total parallel momentum {Pztot) versus total charge (Zfot) for the 
93 Mg system at 30 MeV/nucleon. 

In order to select well-detected events, we used the Pztot — Ztot correlations (see 
figure [1]), the total parallel momentum versus the total charge on an event by 
event basis. For the analysis, only the events with a total parallel momentum 
greater than 60% of the beam momentum were kept, Pztot > 13200 MeV/c. 
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^ i\/f, , ^ 

^ -'"tot ^ 


^ (^ijes ^ till/ Ufa 




Tl 


72.9 ± 10.1 


8.1 ±0.8 


6.6 ±0.2 


34.1 ±2.3 


T2 


120.9 ± 10.6 


11.2 ±0.9 


6.4 ±0.2 


31.9 ±2.7 


T3 


176.4 ± 13.8 


13.4 ± 1.0 


6.3 ±0.2 


30.5 ±2.8 



Table 1 

Average values of the total transverse energy, the total multiplicity, the residue 
velocity and charge for the 3 regions (see text). 

The events were sorted as a function of two observables correlated with the 




impact parameter [30]: E^f, the total transverse energy, and Mtot, the total 
charged particle multiplicity. Three regions are defined according to cuts made 
perpendicular to the ridge line of this correlation, which we refer to as Tl, T2 
and T3 (see figure [2]). The cuts were made in order to retain approximately 
the same number of events in each region. Events in the Tl region correspond 
to less dissipative (peripheral) collisions, while those in T3 correspond to more 
dissipative (central) collisions. Table [1] shows the corresponding average values 
of the total transverse energy, the total multiplicity, the residue (the biggest 
fragment) velocity and charge in the three regions. 
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3.2 Source Reconstruction 



The source (QP or incomplete fusion) is isolated by a selection in parallel 
velocity of different particles, as shown in figure [31 The value of the velocity 
cut was adjusted, for each region and for particles with Z values ranging from 
1 to 10, in order to symmetrize the forward-backward emission in the source 
reference frame. Figure H] shows the result for the region Tl : the energy 
spectra of the particles emitted in the forward or backward hemisphere of the 
source frame are very similar in most cases. Similar features are observed for 
the regions T2 and T3. 






Fig. 3. Invariant velocity diag rcLins Vper versus ^par for protons of tlie regions Tl, 
T2 and T3. 

Table [2] shows the values of the velocity cuts for each particle for each region. 
The charge Zs of the source is obtained as the sum of the charge of all particles 
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Table 2 

Velocity cuts for the source selection for the three regions, 
having a parallel velocity Vpar greater than the cut velocity Vcut- 
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Fig. 4. Backward (black circles)-forward(open circles) energy spectra in the source 
reference frame of the region Tl for different nuclear species ranging from protons 
to fragments {Z > 10). 

Each source velocity component k = x,y, z is evaluated as an average weighted 
by the charge of each particle i of the source : 



' a 



sr^npart T^fc 



(2) 



Since INDRA does not allow a complete mass determination, we made the 
hypothesis that the source has the same isotopic ratio as the projectile : 

As = ^^Zs (3) 
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The neutron multiplicity, nineut, is estimated as the difference between the 
source mass and the sum of the Ai masses of the different ejectiles : 



npart 

rrineut = As - ^ Ai (4) 



where the Ai are estimated with the help of two parameterizations : one from 
a cubic polynomial fit on the /9-valley stabihty : 

Astabie = 1.867Z + 0.016^2 - 1.07 * IQ-^Z^ (5) 



and the other as in references [I6l and [31] : 

Acharity = 2.08Z + 2.9 * lO'^^Z^ (6) 

The excitation energy is obtained by calorimetric recontruction : 



In the equation [7] E^"^"- is the kinetic energy of each particle in the source 
reference frame, Aj (A5 for the source) is the mass defect and -E^e"* the 
neutron kinetic energy proportional to the source temperature : 

E'^elt = 2aT (8) 



where a = 0.75 takes into account the temperature decrease as the source 
emits neutrons [32] . 

In table [3] are reported the average characteristics of the reconstructed source 
for the three regions : the charge Zs, the mass As, the average multiplici- 
ties for protons Mp, alphas M^, fragments Mf^ag (biggest one included); the 
neutron multiplicity M^^^l^ and excitation energy obtained with AgtaUe mass 
parameterization; the neutron multiplicity M^^^^' and the excitation energy 
obtained with Acharity mass parameterization. As it can be seen from the ta- 
ble, the two sets of excitation energies differ by about 0.5 MeV/nucleon. This 
is mainly due to the difference in the deduced neutron multiplicity with the 
two possible choices for the residue mass. Although the average fragment mul- 
tiplicity is low, the excitation energy for the regions T2 and T3 is about ~ 3 
MeV/nucleon, value at which experimental results agree to set the multifrag- 
mentation threshold 
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Table 3 

Characteristics of the reconstructed source : charge, mass, Hght particle, fragment 
and neutron multiplicities and excitation energy. 

4 GEMINI and SMM Backtracing Simulations 

We now compare the data with the results of two statistical codes, GEMINI 
[16] and SMM [17]. One expects that experimental results are relative to an 
ensemble of various sources (defined by distributions of charge, mass and ex- 
citation energy). We therefore applied a backtracing procedure ([I3] -[IS]) to 
GEMINI and to SMM simulations in order to determine these distributions. 

The procedure begins by first choosing the observables which will be used to 
constrain the source distributions. In our case we chose the multiplicity dis- 
tributions of protons, alphas and fragments and the charge distribution of the 
residue and of all the products. This choice was made because these observ- 
ables do not depend on the hypothesis taken for the mass of the source. In a 
second step, the procedure explores a 4-dimensional space constituted by the 
source variable parameters : charge Zq, mass Aq, excitation energy eo and an- 
gular momentum Lq (for GEMINI) or rotational energy Eq°^ (for SMM). The 
source parameter distributions are extracted from this 4-dimensional space 
with the help of the Kolmogorov-Smirnov test applied to simulated and ex- 
perimental observable distributions. The simulated events are filtered in each 
step using a software rephca of the experimental apparatus and conditions. 
This iterative procedure converges to give the parameter distributions of the 
simulated sources, which can then be compared with those reconstructed ex- 
perimentally in order to test the agreement of data with both models. 

The region Tl, representing the most peripheral events, could be affected by 
an experimental bias due to the trigger condition which suppresses low multi- 
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z 

Fig. 5. Comparison between data (full points) and GEMINI results (open squares) 
for region T3. In the top row : Mp : proton multiplicity and Ma : alpha multiplicity; 
in the middle row : Mf : fragment multiplicity and the residue charge distribution. 
In the bottom row : the total charge distribution. 



plicity events. Conversely, the region T3, representing more dissipative events, 
is supposed to be less altered. Therefore we focused our attention on T3, this 
latter being of major interest for the study of multifragmentation threshold. 
Figure [5] shows the comparison between experimental observables (full circles) 
and backtraced observables from GEMINI (open squares) for region T3. The 
charge distribution is poorly reproduced, with a lack of fragment production 
with Z = 6 — 20 together with a strong underestimation of the alpha particle 
multiplicity. In this case the average angular momentum from the GEMINI 
backtracing is about 30 — 35 h, sufficiently high to permit complex particle 
emission. This spin value is in agreement with other findings on rather similar 
systems [33] . 
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The maximum of the source charge distribution, shown in Figure [6l is shifted 
to lower Z values as compared to the experimental data, although both dis- 
tributions have the same mean value (see Table H]). 
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Fig. 6. Source charge distributions : experimental (full points), GEMINI backtraced 
filtered (open squares) for region T3. 




Fig. 7. Source excitation energy for AstaUe (left) and Acharity (right) mass 
parametrization. GEMINI backtraced (open squares) and data excitation energy 
(full points). Region T3. 
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Figure [7] displays the excitation energy distributions for Acharity and AstaUe 
parameterizations. It is possible to note a better agreement for the data treated 
with the Acharity mass parameterization since the mean values of the compared 
distributions are coincident. The comparison with the data according to the 
other mass parameterization shows a shift of about 0.5 MeV/nucleon. 

Figures [HlfTOl present the same comparison as before for the SMM calcula- 
tion. A much better agreement is observed : the residue charge and the total 
charge distributions are nicely reproduced and the disagreement on the alpha 
multiplicity (as well as in the fragment multiplicity) is less severe than with 
GEMINI. 




Fig. 8. Comparison between data (points) and SMM results (squares) for region T3. 
In the top row : Mp : proton multiplicity and Ma '■ alpha multiplicity; in the middle 
row : Mf : fragment multiplicity and the residue charge distribution. In the bottom 
row: the total charge distribution. 

It is worth noting that the best agreement with data was found for SMM 
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calculations performed with a freeze out volume ofV = 2Vo, and a rotational 
energy corresponding to an angular momentum of 35 h for the QP source. 
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Fig. 9. Source charge distributions : experimental (full points) and SMM backtraced 
results (open squares) for region T3. 
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Fig. 10. Source excitation energy for Agtabie (left) and Acharity (right) mass parame- 
terization. SMM backtraced (open squares) and data excitation energy (full points). 
Region T3. 

Figures [9] and [in] show the comparison for the source variables : they display 
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a good agreement between the simulated and the experimental distributions. 
The data treated with the AgtaUe mass parameterization show the same shift 
in the mean value as observed with GEMINI. 







^liable MeV/nucleon 


echariij/MeV/nucleon 


As 


Data 


45. ± 2. 


3.4 ±0.5 


3.8 ±0.4 


101. ±4. 


GEMINI 


44. ± 2. 


3.8 ±0.5 


3.8 ±0.5 


103. ±4. 


SMM 


45. ± 2. 


3.7 ±0.5 


3.7 ±0.5 


104. ±3. 



Table 4 

Data and filtered average values for the source variables obtained from GEMINI 
and SMM backtracing for region T3. 



In table H] are reported the mean values of the source variables for the data, 
for GEMINI and for SMM. Both models give a fair agreement with data. 
One should note that mean values are not sufficient to quantify the agree- 
ment between the experiment and the model. Indeed, the examination of fig- 
ures |5] to [To] shows that overall SMM gives a better description of the data 
than GEMINI. Therefore we can conclude that for excitation energies above 
3 MeV/nucleon, the introduction of some degree of multifragmentation in the 
decay of hot nuclei is mandatory in order to reproduce experimental results. 



5 Discussion and Conclusions 

We also deduced the average source temperatures from both models in order 
to compare them to the current systematics provided by Natowitz et al. in 
reference |llj . 

In table O are shown the average temperature and mass values of the systems 
for each region. A close comparison of these values for the region T3 with those 
obtained for the limiting temperature of reference [11] shows a good agreement 
: the systems of this region may undergo multifragmentation. This is consistent 
with the conclusions drawn in the previous section, i.e. the presence of events 
to be described by a multifragmentation scenario. 





Aexp 


Tcem MeV 


TsMM MeV 


Tl 


91.2 ±9.4 


3.55 




T2 


96.0 ±9.8 


4.54 


4.87 


T3 


101.3 ±9.3 


5.19 


5.56 



Table 5 

Mean temperature and mass values for the sources of the three regions. 
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In perspective one should understand why the alpha particle multiplicity is 
underestimated by GEMINI. This could be connected to the dependence of 
the level density parameter a from the excitation energy ([31],[3S])- In the 
present work, the Lestone parameterization was used , but there are different 
formulae describing this dependence of a that could be used alternatively in 
the backtracing. 

In order to better enhance the differences between the two models, an attempt 
could also be done using the backtracing of dynamical observables , more sen- 
sitive to the Coulomb effects [SB] . One should however recall that some of the 
dynamical observables (like velocities for example) are mass dependent. This 
addresses the problem of the choice of the mass parameterization previously 
evoked. On this subject, we may hope to gain some useful insight by upgrading 
in the mass identification of existing 47r devices like CHIMERAPS or by new 
generation projects like FAZIA (Four 47r A Z Ion Array) [57] . 
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